Exact Study of the Effect of Level Statistics in Ultrasmall Superconducting Grains 
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The reduced BCS model that is commonly used for ultrasmall superconducting grains has an 
exact solution worked out long ago by Richardson in the context of nuclear physics. We use it to 
check the quality of previous treatments of this model, and to investigate the effect of level statistics 
on pairing correlations. We find that the ground state energies are on average somewhat lower for 
systems with non-uniform than uniform level spacings, but both have an equally smooth crossover 
from the bulk to the few-electron regime. In the latter, statistical fluctuations in ground state 
energies strongly depend on the grain's electron number parity. 



PACS number: 74.20. Fg, 74. 25. Ha, 74.80. Fp 

Recent experiments by Ralph, Black and Tinkham, in- 
volving the observation of a spectroscopic gap indicative 
of pairing correlations in ultrasmall Al grains have 
inspired a number of theoretical f^|li"|] studies of how su- 
perconducting pairing correlations in such grains are af- 
fected by reducing the grains' size, or equivalently by in- 
creasing its mean level spacing d oc VoP 1 until it exceeds 
the bulk gap A. In the earliest of these, a grand-canonical 
(g.c.) BCS approach was applied to a reduced BCS 
Hamiltonian for uniformly spaced, spin-degenerate lev- 
els; it suggested that pairing correlations, as measured 
by the condensation energy E c , vanish abruptly once d 
exceeds a critical level spacing d c that depends on the 
parity (0 or 1) of the number of electrons on the grain, 
being smaller for odd grains (d\ ~ 0.89A) than even 
grains (<i§ ~ 3.6A). A series of more sophisticated canon- 
ical approaches (summarized below) confirmed the par- 
ity dependence of pairing correlations, but established 
that the abrupt vanishing of pairing correlations 
at d c is an artifact of g.c. treatments: pairing correla- 
tions do persist, in the form of so-called fluctuations, 
to arbitrarily large level spacings, and the crossover be- 
tween the bulk superconducting (SC) regime (d <C A) 
and the fluctuation-dominated (FD) regime (d 3> A) is 
completely smooth jl(J. Nevertheless, these two regimes 
are qualitatively very different the condensation 

energy, e.g., is an extensive function of volume in the for- 
mer and almost intensive in the latter, and pairing cor- 
relations are quite strongly localized around the Fermi 
energy ep, or more spread out in energy, respectively. 

After the appearance of all these works, we became 
aware that the reduced BCS Hamiltonian on which they 
are based actually has an exact solution. It was published 
by Richardson in the context of nuclear physics (where it 
is known as the "picket-fence model" ) , in a series of pa- 
pers between 1963 and 1977 |1^,[f|] which seem to have 
completely escaped the attention of the condensed mat- 
ter community. The beauty of this solution, besides its 
mathematical elegance, is that it also works for the case 



of randomly-spaced levels. It thus presents us with two 
rare opportunities, which are the subject of this Letter: 
(i) to compare the results of various previously-used ap- 
proximations against the benchmark set by the exact so- 
lution, in order to gauge their reliability for related prob- 
lems for which no exact solutions exist; and very inter- 
estingly, (ii) to study the interplay of randomness and 
interactions in a non-trivial model exactly, by examining 
the effect of level statistics on the SC/FD crossover. 

There is a previous study of the latter question by 
Smith and Ambegaokar using the g.c. mean-field BCS 
approach (^], who concluded, interestingly, that random- 
ness enhances pairing correlations: compared to the case 
of uniform spacings |2j , they found that a random spacing 
of levels (distributed according to the gaussian orthogo- 
nal ensemble) on average lowers the condensation energy 
E c to more negative values and increases the critical level 
spacings at which E c vanishes abruptly, but these still 
are parity dependent ((dj) = 1.8A, (d§) ~ 14A). How- 
ever, the abrupt vanishing of E c found by SA can be 
suspected to be an artifact of their g.c. mean-field treat- 
ment, as was the case in Indeed, our exact results 
for random levels show (1) that the SC/FD crossover 
is as smooth as for the case of uniformly-spaced levels; 
this means, remarkably, that (2) even in the presence of 
randomness pairing correlations never vanish, no mat- 
ter how large d/A becomes; quite the opposite, (3) the 
randomness-induced lowering of E c is strongest in the 
FD regime; in the latter, moreover, (4) the statistical 
fluctuations in E depend quite strongly on parity. 

Exact solution. — Ultrasmall superconducting grains 
are commonly described PHTl[ by a reduced BCS model, 
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for a set of pairs of time-reversed states \j, ±) with en- 
ergies £j, mean level spacing d and dimensionless cou- 
pling constant A. Unbeknownst to the authors of P~pT|, 
Richardson had long ago solved this model exactly, for 
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an arbitrary set of levels £j (not necessarily all distinct): 
Since singly- occupied levels do not participate in and re- 
main "blocked" to the pairscattering described by 
H , the labels of such levels are good quantum numbers. 
Let \n, B) thus denote an eigenstate with N = 2n + b 
electrons, b of which sit in a set B of singly-occupied, 
blocked levels, thereby contributing £b — J2ieB £i *° * ne 
total energy. The dynamics of the remaining n pairs is 
then governed by 
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where the pair operators bj = Cj—Cj+ satisfy "hard-core 
boson" commutation relations, = 5jf(l — 2bjbj), 

and the sums are over all unblocked levels. Richardson 
showed that the lowest-lying of the eigenstates |n, B) has 
an (unnormalizcd) product form, 
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where the n parameters e v (v — l,...,n) are that par- 
ticular solution of the n coupled algebraic equations 
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that yields the lowest value for the "pair energy" £(n) = 
£™=l e v- Moreover, \n, B)c has total energy £(n) + £b- 
The lowest-lying of all eigenstates with n pairs and b 
blocked levels, say \n, o)g with energy £^{n), is that 
\n,B)c for which the blocked levels in B are all as close 
as possible to sp, the Fermi energy of the uncorrelated 
TV-electron Fermi sea \Fn). In this Letter we shall always 
take all the Sj to be non-degenerate. The e„, which may 
be thought of as self-consistently-determined pair ener- 
gies, then coincide at A = with the lowest n energies 
2ej (j = 1, . . . , n), and smoothly evolve toward lower val- 
ues as A is turned on. This fact can be exploited during 
the numerical solution of @ , which can be simplified by 
first making some algebraic transformations, discussed in 
detail in |13| , that render the equations less singular. 

Uniformly-spaced levels. — Our first application of the 
exact solution is to check the quality of results previ- 
ously obtained by various other methods. Most previ- 
ous works [§-^,^-0 studied a half-filled band with fixed 
width 2uji) of uniformly-spaced levels (i.e. e.j — j d), con- 
taining N — 2n + b electrons. Then the level spacing 
is d = 2uo£,/N and in the limit d — > the bulk gap is 
A = w sinh(l/A) _1 . Following [§, we take A = 0.224 
throughout this paper. To study the SC/FD crossover, 
two types of quantities were typically calculated as func- 
tions of increasing d/A, which mimics decreasing grain 
size: the even and odd (b = 0, 1) condensation energies 
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FIG. 1. (a) The even and odd (b = 0, 1) condensation en- 
ergies E% of Eq. @, calculated with BCS, PBCS and exact 
wave functions, as functions of d/A — 2sinh(l/A)/(2n + b), 
for A = 0.224. For comparison the dotted line gives the 
"bulk" result E$ ulk = -A 2 /(2d). (b) Comparison of the 
parity parameters A ML Q of Eq. (||) obtained by various 
authors: ML's analytical result (dotted lines) [A(l — d/2A) 
for d < A, and d/21og(od/A) for d > A, with a = 1.35 
adjusted to give asymptotic agreement with the exact result]; 
grand-canonical BCS approach (dash-dotted line) [the naive 
perturbative result |Ad is continued to the origin]; PBCS ap- 
proach (short-dashed line); Richardson's exact solution (thick 
solid line); exact diagonalization and scaling by MFF (open 
circles) and BH (long-dashed line). 

and a parity parameter introduced by Matveev and 
Larkin (ML) || to characterize the even-odd ground 
state energy difference, 



A (n) = £?(n) - [#(n) + £^{n + l)]/2 . 



(5) 
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Following the initial g.c. studies [Ert|, the first canonical 
study was that of Mastellone, Falci and Fazio (MFF) 0, 
who used Lanczos exact diagonalization (with n < 12) 
and a scaling argument to probe the crossover regime. 
Berger and Halperin (BH) M showed that essentially the 
same results could be achieved with n < 6 by first re- 
ducing the bandwidth and renormalizing A, thus signifi- 
cantly reducing the calculational effort involved. To ac- 
cess larger systems and fully recover the bulk limit, fixed- 
n projected variational BCS wavefunctions (PBCS) were 
used in (for n < 600); significant improvements over 
the latter results, in particular in the crossover regime, 
were subsequently achieved in JT(| using the density ma- 
trix renormalization group (DMRG) (with n < 400). 
Finally, Dukelsky and Schuck |ll[ showed that a self- 
consistent RPA approach, that in principle can be ex- 
tended to finite temperatures, describes the f.d. regime 
rather well (though not as well as the DMRG). 
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To check the quality of the above methods, we |T^] 
computed (n) and A ML (n) using Richardson's solu- 
tion (Fig. 0). The exact results (a) quantitatively agree, 
for d— > 0, with the leading —A 2 /2d behavior for (n) 
obtained in the g.c. BCS approach |§-f|], which in this 
sense is exact in the bulk limit, corrections being of order 
d°; (b) confirm that a completely smooth |l(]] crossover 
occurs around the scale d ~ A at which the g.c. BCS ap- 
proach breaks down; (c) show that the PBCS crossover 
H is qualitatively correct, but not quantitatively, be- 
ing somewhat too abrupt; (d) are reproduced remarkably 
well by the approaches of MFF Hand BH gj; (e) are 
fully reproduced by the DMRG of with a relative er- 
ror of < 10~ 4 for n < 400; our figures don't show DMRG 
curves, since they are indistinghuishable from the exact 
ones and are discussed in detail in fl(i|| . 

The main conclusion we can draw from these compar- 
isons is that the two approaches based on renormalization 
group ideas work very well: the DMRG is essentially ex- 
act for this model, but the band-width rescaling method 
of BH also gives remarkably (though not quite as) good 
results with rather less effort. In contrast, the PBCS 
approach is rather unreliable in the crossover region. 

Randomly-spaced levels. — The remainder of this Letter 
addresses the question of how randomness of the levels £j 
affects pairing correlations. We studied half- filled bands 
of ./V = 2n + b non- uniformly spaced but non-degenerate 
levels (for N < 260), with b = 0, 1. The energy levels in 
small metallic grains with time reversal symmetry follow 
the gaussian orthogonal ensemble distribution [jL6| . We 
generated sets of levels e j (i = 1 , • • • , N) by diagonalizing 
2Nx2N random matrices, taking N adjacent values from 
the central part of the eigenspectrum (to avoid boundary 
effects) and performing the rescaling 
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to ensure an average level spacing of one in units of d. In 
Figure |^ we show four such sets of randomly generated 
levels for N = 28, together with the equally spaced set. 

For each such set of 2n + b levels, we calculated the ex- 
act ground state energy £^(n), the condensation energy 
Efr(n), and the spectroscopic gap 0] 
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which gives the energies needed to break a single pair 
in the (even or odd) ground state. Subsequently we 
calculated the ensemble average (E^(n)) and variance 

SE^(n) = ^({E^) 2 ) - (EC) 2 (and analogously (J3f) 

and 5E^ ) over many realizations of random matrices. 
The ensemble size was 1000 for 24 < N < 40, and varied 
between 700 and 150 for 40 < N < 260. Figure | presents 
our results for these ensemble-averages (solid lines, with 
variances indicated by fluctuation bars) together with 



those for the uniformly-spaced (u.s.) set discussed above 
(dashed lines). It shows a number of interesting features. 

Firstly, the two main conclusions of SA [|| are con- 
firmed, namely (a) that pairing correlations are on av- 
erage stronger for randomly- than for uniformly-spaced 
levels, (Eff) < Efr(u.s.); and (b) that the parity effect 
persits in the presence of randomness, (Eq) < (E^). In 
SA's g.c. calculation these facts could be understood || 
from a condition, derived from the BCS gap equation, 
for having non- vanishing pairing correlations, namely 
2/A < J2j^B 1 / 1 ^ j" — Ml- Here Sj and the g.c. chemical 
potential \x are in units of d, and the number of terms in 
the sum is of order 2cur)/d. As d increases, this number 
decreases, until the inequality ceases to hold at a criti- 
cal spacing d%. Since statistical fluctuations to smaller 
values of \e~j — p\ carry more weight than those to larger 
values, fluctuations on average tend to increase which 
explains (a); moreover, since the blocking of levels close 
to /2 reduces the number of terms in the sum, it reduces 
dfr, which explains (b). 

Since the equation on which SA's elegant argument 
is based breaks down in the FD regime, let us attempt 
another way of interpreting (a) and (b): pairing corre- 
lations involve a non-zero amplitude to find pair states 
with £j > £p doubly occupied and ones with e.j < ef 
empty. Such correlations between states below and above 
£f, called "pair- mixing across £f" in @, gain interaction 
energy but cost some kinetic energy. The latter cost is 
the smaller, the closer the states involved in pair-mixing 
across e f he together (which is why the bulk limit d — > 
is so strongly correlated). Statistical fluctuations in level 
positions that yield more-closely or less-closely spaced 
levels around Sf than for the uniform case, would thus 
cause a respectively lower or higher kinetic energy cost 
for pairmixing across sp; according to (a), the former on 
average outweighs the latter, just as had SA concluded 
in H|. Furthermore, in odd grains the blocked level at 
Ef always causes the spacing between pair levels below 
and above ep, and hence the kinetic energy cost for pair- 
mixing across Bp, to be somewhat larger than in even 
grains, which explains (b). 

Now, the ability of the exact solution to correctly 
treat the FD regime enables us to uncover several fur- 
ther facts that are beyond the reach of SA's g.c. mean 
field approach: (c) The SC/FD crossover is as smooth 
for randomly- as for uniformly-spaced levels, confirming 
that the abrupt vanishing of pairing correlations at some 
critical level spacing found by SA is an artifact of their 
g.c. mean field treatment, just as in ^-^|. (d) Even in the 
presence of randomness, pairing correlations never van- 
ish, no matter how large d/ A. Quite the opposite, (e) the 
randomness-induced lowering in condensation energy to 
more negative values, (E^ ) — E^(u.s.), is strongest in 
the FD regime; this perhaps somewhat counterintuitive 
result illustrates that the smaller the number of levels is 
that lie "close to" (i.e. within A of) ep, the stronger is 
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the effect of fluctuations in their positions on the kinetic 
energy cost for pair-mixing; conversely, this randomness- 
induced lowering of E^ decreases in the crossover regime 
and becomes negligible in the SC regime, in which very 
many levels lie within A of Ep- (f) The variances SE^ 
are essentially (^-independent in the range 24 < N < 260, 
implying that the relative statistical fluctuations of E^ 
should be negligible in the bulk limit, as expected. 

Remarkably, we can also discern (g) three "parity- 
dependent fluctuation effects" , in that the following three 
quantities are larger for even than for odd grains: (gl) 
the variances SE^ (with SEfi ~ 25E? ~ A/2); and 
the randomness-induced changes in (g2) condensation 
energies \(E£) — E^'(u.c)\ and (g3) spectroscopic gaps 
\(Eg ) - £ h G (u.s.)| (inset of Fig. g). All three of these 
effects have the same origin as the more familiar par- 
ity effect (b), namely blocking: the more levels around 
Ef are blocked, the larger the effective spacing between 
states involved in pair- mixing across £f, and hence the 
smaller the sensitivity of the total energy to statistical 
fluctuations in level positions. 

In conclusion, using Richardson's exact solution we 
have found that level randomness does not modify the 
smooth nature of the SC/FD crossover. It just en- 
hances pairing correlations somewhat compared to those 
of uniformly-spaced levels, having the strongest effect in 
the FD regime. In the latter we found that statistical 
fluctuations become strongly parity dependent. 
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FIG. 2. Sets of energy levels with TV 
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28. Set c has equally 



spaced levels, with spectroscopic gap [Eq. (|7j)] E§ /d = 1.54. 
Sets a, b (or c, d) are randomly spaced; among all sets with 
N — 28 we studied, the ones shown have the smallest (largest) 
values for E$ /d, namely 0.886, 0.891 (3.30,3.37), due to the 
small (large) spacing between the two levels closest to ef, 
illustrating how random level fluctuations affect energy gaps. 
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FIG. 3. Exact even and odd condensation energies, E% for 
equally spaced levels (dashed line), and the ensemble-average 
(Et) for randomly-spaced levels (solid line). The height of 
the fluctuation bars gives the variances SE^ . The inset shows 
the corresponding spectroscopic gaps E® and variances SE^ . 
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